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The evolution of the magnetic moment in iron clusters containing 20 to 400 atoms is investigated 
using first-principles numerical calculations based on density-functional theory and real-space pseu- 
dopotentials. Three families of clusters are studied, characterized by the arrangement of atoms: 
icosahedral, body-centered cubic centered on an atom site, and body-centered cubic centered on 
the bridge between two neighboring atoms. We find an overall decrease of magnetic moment as the 
clusters grow in size towards the bulk limit. Clusters with faceted surfaces are predicted to have 
magnetic moment lower than other clusters with similar size. As a result, the magnetic moment is 
observed to decrease as function of size in a non-monotonic manner, which explains measurements 
performed at low temperatures. 



The existence of spontaneous magnetization in metal- 
lic systems is an intriguing problem because of the ex- 
tensive technological applications of magnetic phenom- 
ena and an incomplete theory of its fundamental mech- 
anisms. Clusters of metallic atoms are important in this 
respect as they serve as a bridge between the atomic limit 
and the bulk, and can form a basis for understanding the 
emergence of magnetization as function of size. Several 
phenomena such as ferromagnetism, metallic behavior, 
and ferroelectricity have been intensely explored in bulk 
metals, but the way they manifest themselves in clusters 
is an open topic of debate. At the atomic level, ferro- 
magnetism is associated with partially filled 3d orbitals. 
In solids, ferromagnetism may be understood in terms 
of the itinerant electron model Q, which assumes par- 
tial derealization of the 3d orbitals. In clusters of iron 
atoms, derealization is weaker owing to the presence of 
a surface, whose shape affects the magnetic properties 
of the cluster. Because of their small size, iron clusters 
containing a few tens to hundreds of atoms are superpara- 
magnetic: the entire cluster serves as a single magnetic 
domain, with no internal grain boundaries 2]. Conse- 
quently, these clusters have strong magnetic moments, 
but exhibit no hysteresis. 

The magnetic moment of nano-sized clusters has been 
measured as function of temperature and size 0, 0, Q , 
and several aspects of the experiment have not been 
fully clarified, despite the intense work on the subject 
0- Q- S H \Hi El] • One intriguing experimental observa- 
tion is that the specific heat of such clusters is lower than 
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the Dulong-Petit value, which may be due to a magnetic 
phase transition p|. In addition, the magnetic moment 
per atom does not decay monotonically as function of 
the number of atoms and for fixed temperature. Possi- 
ble explanations for this behavior are: structural phase 
transitions, strong dependence of magnetization with the 
shape of the cluster, or coupling with vibrational modes 
|4|. One difficulty is that the structure of such clusters 
is not well known. First-principles and model calcula- 
tions have shown that clusters with up to 10 or 20 atoms 
assume a variety of exotic shapes in their lowest-energy 
configuration 0,0. For larger clusters, there is evi- 
dence for a stable body-centered cubic (BCC) structure, 
which is identical to ferromagnetic bulk iron p. 

The evolution of magnetic moment as function of clus- 
ter size has attracted considerable attention [S 0, IS IS 
S IS E3 ES El El El A key question to be resolved 
is: What drives the suppression of magnetic moment as 
clusters grow in size? In the iron atom, the permanent 
magnetic moment arises from exchange splitting: the 3d^ 
orbitals (majority spin) are lower in energy and com- 
pletely occupied with 5 electrons, while the 3d^ orbitals 
(minority spin) are partially occupied with one electron, 
resulting in a magnetic moment of 4/i#, \±b being the 
Bohr magneton. When atoms are assembled in a crystal, 
atomic orbitals hybridize and form energy bands: 4s or- 
bitals create a wide band which remains partially filled, 
in contrast with the completely filled 4s orbital in the 
atom; while the 3d^ and 3d^ orbitals create narrower 
bands. Orbital hybridization together with the differ- 
ent bandwidths of the various 3d and As bands result in 
weaker magnetization, equivalent to 2.2 /i^/atom in bulk 
iron. 

In atomic clusters, orbital hybridization is not as 
strong because atoms on the surface of the cluster have 
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FIG. 1: Density of states in the clusters Fe38s(a) and Fe393(b), 
majority spin (upper panel) and minority spin(lower panel). 
Fe388 corresponds to a fragment of the BCC crystal. Fe393 
has local icosahedral coordination. For reference, the density 
of states in bulk iron is shown in dashed lines. The Fermi 
energy is chosen as energy reference. 



fewer neighbors. The strength of hybridization can be 
quantified by the effective coordination number. A the- 
oretical analysis of magnetization in clusters and thin 
slabs indicates that the dependence of the magnetic mo- 
ment with the effective coordination number is approx- 
imately linear || 0, . But the suppression of mag- 
netic moment from orbital hybridization is not isotropic 
[l4| . For instance, if we consider a layer of atoms on the 
x — y plane, the 3d orbitals oriented on the plane (the 
ones with angular dependence xy and x 2 — y 2 ) will hy- 
bridize more effectively than orbitals oriented normal to 
the plane (with angular dependence 3z 2 — r 2 ), because 
of increased overlap among orbitals on the x — y plane. 
If the layer has infinite extent, orbitals oriented on that 
plane will tend to form an energy band wider than the 
one formed by orbitals normal to it. With a wider band, 
the difference in electron population between the major- 
ity spin channel and the minority spin channel is reduced, 
leading to reduced spin polarization and overall weaker 
magnetization. Reduced orbital hybridization has been 
observed to enhance the magnetic moment of metallic 
clusters and thin layers [H| . As a consequence of this 
anisotropy, clusters with faceted surfaces are expected to 
have magnetic properties different from clusters with ir- 
regular surfaces, even if they have the same effective co- 



ordination number. This effect is likely responsible for a 
non-monotonic suppression of magnetic moment as func- 
tion of cluster size. In order to analyze the role of surface 
faceting, we have performed first-principles calculations 
of the magnetic moment of iron clusters with various ge- 
ometries and with sizes ranging from 20 to 400 atoms. 

We determine the electronic structure of clusters 
within the framework of pseudoptentials|l7| constructed 
using density-functional theory (DFT) [IlElIi^- DFT 
is an established theory for first-principles studies of 
weakly and moderately correlated electronic systems. 
The exchange-correlation functional used in this work 
employs the generalized gradient approximation (GGA) 
|2lj . We have observed that the GGA predicts magnetic 
moments enhanced with respect to the simpler local- 
density approximation by 2 % to 10 %. For a fixed geom- 
etry of the cluster, we solve self-consistently the Kohn- 
Sham equation on a regular grid in real space Hj| . 
Proper boundary conditions are obeyed by imposing the 
electronic wave functions to vanish on the boundary of 
a large spherical domain, which contains the system of 
interest. No explicit basis set is used. Numerical con- 
vergence is controlled with two parameters: the radius of 
the domain (typically 5 A larger than the cluster radius) 
and the nearest-neighbor spacing in the regular grid. For 
iron atoms, we use a spacing of 0.3 a.u., (approximately 
0.16 A). 

We use the PARSEC codejUH^. This code makes use 
of symmetry properties of the system and very efficient 
techniques for solving the Kohn-Sham equation^? E3- 
Until recently, a significant fraction of numerical effort 
was spent in perform ing exact diagonalization of the 
Kohn-Sham equation [2Q, E2, EJ EJ- Curently, this 
step is replaced with a series of subspace filtering iter- 
ations with Chebyshev polynomials, which reduce the 
overall numerical effort by one order of magnitude or 
more|24|. This dramatic advance in methodology allows 
us to study confined systems with hundreds, if not thou- 
sands of atoms in a straightforward and computationally 
efficient manner. 

Clusters of both icosahedral and BCC symmetry are 
explored in our work. In order to investigate the role of 
surface faceting, we construct clusters with faceted and 
non-faceted surfaces. Faceted clusters are constructed 
by adding successive atomic layers around a nucleation 
point. Small faceted icosahedral clusters exist with sizes 
13, 55, 147, and 309. Faceted BCC clusters are con- 
structed with BCC local coordination and, differently 
from icosahedral ones, they do not need to be centered 
on an atom site. We consider two families of cubic clus- 
ters: atom centered or bridge centered, respectively for 
clusters with nucleation point at an atom site or on the 
bridge between two neighboring atoms. The lattice pa- 
rameter is equal to the bulk value, 2.87 A. Non-faceted 
clusters are built by adding shells of atoms around a nu- 
cleation point so that their distance to the nucleation 
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FIG. 2: Magnetic moment versus spin moment calculated for 
the atom-centered BCC ("plus" signs) and bridge-centered 
BCC (crosses) iron clusters. The approximate ratio is 
M/(S z )=g e ff = 2MfjL B /h. 



point is less than a specified value. As a result, non- 
faceted clusters usually have narrow steps over otherwise 
planar surfaces and the overall shape is almost spherical. 
By construction, non-faceted clusters have well-defined 
point-group symmetries: Ih or Th for the icosahedral 
family; Oh for the atom-centered family; and for 
the bridge-centered family. Clusters constructed in that 
manner show low tension on the surface, making surface 
reconstruction less likely. Our calculations indicate that 
atoms on the surface feel forces weak in magnitude and 
directed towards the center of the cluster. Owing to the 
small surface/ volume ratio in large clusters, the impact 
of surface reconstruction on those clusters will be small, 
if not negligible. 

As clusters grow in size, their properties approach the 
properties of bulk iron. Figure da) shows the density 
of states (DOS) for Fe388, with local BCC coordination. 
At this size range, the density of states assumes a shape 
typical of bulk iron, with a three-fold partition of the 3d 
bands. In addition, the cohesive energy of this cluster is 
only 77 meV lower than in bulk. This evidence suggests 
that interesting size effects will be predominantly ob- 
served in clusters smaller than Fe388- Figure [^b) shows 
the DOS for Fe393, which belongs to the icosahedral fam- 
ily. This cluster has a very smooth DOS, with not much 
structure compared to Fe388 and bulk BCC iron. This 
is due to the icosahedral-like arrangement of atoms in 
Fe393. The overall dispersion of the 3d peak ( 4 eV for 
3d^ and 6 eV for 3d;) is nevertheless similar in all the 
calculated DOS. 

The magnetic moment is calculated as the expectation 
value of the total angular momentum: 



M = [ f[g s (S z ) 
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FIG. 3: Calculated magnetic moments for clusters in the 
atom-centered ("plus" signs, a), bridge-centered (crosses, b), 
and icosahedral (triangles, c) families. Experimental data |3 
is shown in black diamonds with error bars. Some of the 
faceted and non-faceted clusters are depicted next to their cor- 
responding data points. The dashed lines indicate the value 
of magnetic moment per atom in bulk iron. 



where g s = 2 is the electron gyromagnetic ratio. Figure [2| 
illustrates the approximately linear dependence between 
the magnetic moment and spin moment, (S z ), through- 
out the whole size range. This results in an effective 
gyromagnetic ratio g e ff = 2.04/ie/Ti, which is somewhat 
smaller than the gyromagnetic ratio in bulk BCC iron, 
2.09/i#/ft. The difference in ratios is probably due to 
an underestimation in the orbital contribution, (L z ). In 
the absence of an external magnetic field, orbital mag- 
netization arises from the spin-orbit interaction, which is 
included in the theory as a model potential, 



V*< 



-£L • S 



(2) 



(1) 



where £ = 80 meV/ft 2 0. 

Figure|3|shows the magnetic moment of several clusters 
belonging to the three families studied: atom-centered 
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BCC (top panel), bridge- centered BCC (middle panel) 
and icosahedral (bottom panel). Experimental data ob- 
tained by Billas and collaborators 3] is also shown. The 
suppression of magnetic moment as function of size is 
readily observed. Also, clusters with faceted surfaces are 
predicted to have magnetic moments lower than other 
clusters with similar sizes. This is attributed to more ef- 
fective hybridization of d orbitals along the plane of the 
facets. We also notice that the correlation between mag- 
netic moment and surface smoothness is not always well 
defined. For instance, the cluster Fei75 m FigureOJa) has 
wide facets but its magnetic moment is not much weaker 
than clusters with similar sizes. In Figure Efc), clusters 
with strongly suppressed magnetic moment have 67 and 
135 atoms, whereas clusters with wide facets and similar 
sizes have 55 and 147 atoms. 

The measured non-monotonic behavior of magnetic 
moment can be attributed to the shape of the surface. 
Under this assumption, islands of low magnetic moment 
(observed at sizes 45, 85 and 188) are associated to clus- 
ters with faceted surfaces. In the icosahedral family, 
the islands of low magnetic moment are located around 
faceted clusters containing 55, 147, and 309 atoms. The 
first island is displaced by 10 units from the measured lo- 
cation. For the atom-centered and bridge-centered fam- 
ilies, we found islands at (65, 175) and (92, 173) respec- 
tively, as indicated in FigureOJa,b). The first two islands 
are also close to the measured islands at 85 and 188. 
Clearly, there is no exact superposition in the location of 
calculated islands and measured islands. The magnetic 
moment was measured in clusters at 120 K 0, U- At 
that temperature, vibrational modes or the occurrence 
of metastable configurations can shift the islands of low 
magnetic moment or make them more diffuse. Assum- 
ing that the non-monotonic decay of magnetic moment 
is dictated by the cluster shape, we also conclude that 
clusters with local structures different from the ones we 
discuss here (such as cobalt clusters with hexagonal-close 
packed coordination, or nickel clusters with face-centered 
cubic coordination) should have islands of low magnetic 
moment located at different "magic numbers" , according 
to the local atomic coordination. 

In summary, we discussed the behavior of magnetiza- 
tion in iron clusters containing 20 to 400 atoms in the 
light of first-principles density-functional theory. The 
magnetic moment is found to decay as function of clus- 
ter size in a non-monotonic fashion: clusters with faceted 
surfaces are predicted to have magnetic moments lower 
than clusters with non-faceted surfaces. As a conse- 
quence, clusters with many steps, or atoms protruding 
from the surface, are expected to have strong magnetic 
properties at low temperatures. In addition, large clus- 
ters with icosahedral structure are expected to have mag- 
netic moments lower than clusters with BCC structure. 
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